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Abstract. High-luminosity e + e~ colliders at the GeV scale (flavor factories) have been recently recognized 
to be an ideal environment to search for a light weakly coupled vector boson U (dark photon) emerging 
in several new physics models. At flavor factories a particularly clean channel is the production of the U 
boson in association with a photon, followed by the decay of the U boson into lepton pairs. Beyond the 
approximations addressed in previous works, we perform an exact lowest order calculation of the signal and 
background processes of this channel. We also include the effect of initial and final state QED corrections 
neglected so far, to show how they affect the distributions of experimental interest. We present new results 
for the expected statistical significance to a dark photon signal at KLOE/KLOE-2 and future super- B 
factories. The calculation is implemented in a new release of the event generator BabaYaga@NLO, which 
is available for full event simulations and data analysis. 

PACS. 12.60.Cn Extensions of electroweak gauge sector - 14.70.Pw Other gauge bosons - 13.66.De Lepton 
production in e~e + interactions - 13.40.Ks Electromagnetic corrections to strong- and weak-interaction 
processes 
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1 Introduction 

In recent years, striking astrophysical observations have 
failed to find an interpretation in terms of standard as- 
trophysical or particle physics sources. Among these ob- 
servations, there are the 511 KeV gamma-ray signal from 
the galactic center observed by the INTEGRAL satellite 
[1 , the excess in the cosmic ray positrons reported by 
PAMELA [2], the total electron and positron flux mea- 
sured by ATIC [3J, Fermi [4 , and HESS [5], the annual 
modulation of the DAMA/LIBRA signal [6] and the fea- 
tures of the low-energy spectrum of rare events reported 
by the CoGeNT collaboration [3[8]. These evidences can 
be comprehensively interpreted by well motivated exten- 
sions of the Standard Model (SM) that predict the exis- 
tence of a WIMP Dark Matter (DM) particle belonging to 
a secluded gauge sector under which the SM particles are 
uncharged [9lfT( H [TT j ri 2 [1 1 3 [ HH] . In particular, an abelian 
gauge symmetry, with an associated U boson ( "dark pho- 
ton"), can communicate with the SM through a kinetic 
mixing term of the form 



r ■ — 

^mix — 



" 9 \iv ^dark 
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where e is a kinetic mixing parameter. Annihilation of 
DM into the U boson, which decays into charged lep- 
tons and is light enough to kinematically forbid a decay 



that produces antiprotons, can explain the electron and/or 
positron excesses, and the absence of a similar effect in 
the PAMELA antiproton observations. Independently of 
connections with DM physics, general theoretical argu- 
ments P^fTB] suggest that the kinetic mixing parameter 
e must be naturally of the order of 10 -4 — 10 -2 , if the 
mixing occurs through loop effects. On the other hand 
the phenomenological condition Mjj below the GeV scale, 
driven by the astrophysical observations described above, 
requires that e is in the same range in order to avoid con- 
tradiction with the available data, in particular with the 
constraints imposed by the precision measurements of the 
anomalous magnetic moment of the leptons and of the 
electromagnetic coupling constant [T7] . 

From a phenomenological perspective, the motivations 
just discussed are, perhaps, the most compelling towards 
postulating the existence of a new gauge boson U. It is 
however worth mentioning that there is a broad class of 
new physics models that predict the existence of such 
a particle, without relying on DM arguments (see e.g. 
Refs. HS1IIS1). 

An intriguing consequence of the above ideas is that 
such a light U boson (if it exists) can be directly pro- 
duced in a controlled environment, such as fixed target 
experiments [18,19,20 or high-luminosity e + e~ colliders 
at the GeV scale (flavor factories) [2ll[22j[23l[2l[25l[2i[271 
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128] . The status and perspectives of these experiments are 
the topics studied in the SLAC workshop "Searches for 
New Forces at the GeV-scale" [29]. At flavor factories, e.g. 
at DA^NE, BESIII and present and future ^-factories, a 
particularly clean and simple channel, which is insensitive 
to the details of the Higgs sector of the secluded group, 
is the associated production of a U boson and a photon, 
with decay of the U into lepton pairs [21,23,26,28 . A dis- 
tinctive feature of the expected signal is the appearance 
of a Breit-Wigner peak in the shape of the invariant mass 
distribution of the lepton pairs induced by the mechanism 
of photon radiative return and corresponding to U bo- 
son resonant production. The drawback of this channel is 
the fairly small value, over a wide range of the param- 
eters, of the signal cross section in comparison with the 
rate of the large physics backgrounds given by the QED 
radiative processes e + e _ — »■ / + /~7, I = e,/i, which can 
be rejected by cutting on the invariant mass of the lepton 
pair. The studies present in the literature agree on the con- 
clusion that the Uj production process allows to reach a 
sensitivity to the kinetic mixing parameter in the range 
e ~ 10 -3 — 10 -2 at present flavor factories, for a U boson 
with mass M\j up to a few GeV. These analyses are gener- 
ally based on the evaluation of the number of signal events 
through the calculation of the differential cross section of 
the 2^2 process e + e~ — )> Uj : including the decay of the 
on-shell U boson into lepton pairs by means of branch- 
ing ratios [21,26,28 , and/or an approximate estimate of 
the backgrounds. More importantly, all the studies so far 
performed neglect the contribution of higher order initial 
and final state QED corrections, which are known to be 
a phenomenologically relevant effect at GeV-scale e + e~ 
colliders [30] . 

The aim of the paper is detailed in the following. First, 
we perform an exact tree level calculation of the signal 
and background processes contributing to the signatures 
e + e~ — »> /i + /i~7, e + e _ 7. Second, we compute the effect of 
the most important higher order corrections induced by 
multiple photon radiation and vacuum polarization, that 
are a source of relevant systematic effects. Third, we use 
our calculation in order to assess its impact on the experi- 
mental sensitivity as evaluated in the literature, and show 
how this can be enhanced by means of event selections 
not considered so far. Last but not least, the calculation 
is made available in the form of a Monte Carlo genera- 
tor, which is still missing for such studies [24,29 and can 
be used for data analysis at flavor factories to study any 
physical observable, i.e. invariant mass distributions but 
also angular distributions, correlations and so on. 

The paper is organized as follows. In Section [2] we 
briefly review the model under consideration and describe 
the details of our calculation. We also present a first sam- 
ple of numerical results, showing in particular how the 
radiative corrections significantly distort the shape of the 
invariant mass distribution which is the most important 
observable for experimental searches. In Section [3] we give 
new results for the statistical significance to a dark pho- 
ton signal at KLOE/KLOE-2 and a super- £? factory, to 
complete the picture existing in the literature for the ex- 




Fig. 1. Examples of Feynman diagrams with dark photon ex- 
change contributing to the process e + e~ — » 7, U — ► / + /~7, 
/ = e, /x. The total number of diagrams is eight for the fi final 
state, sixteen for the e final state, where also t-channel contri- 
butions exist. 



perimental sensitivity of flavor factories. In Section H] we 
draw our conclusions. 



2 Theoretical Framework 

Following most of the phenomenological works in the lit- 
erature, we consider a minimal implementation of a se- 
cluded 11(1)5 sector [22]. The SM Lagrangian is modi- 
fied by the inclusion of a 11(1)5 gauge group which con- 
tains DM fields, a vector gauge field A'^ and a single com- 
plex scalar Higgs field <\> responsible for spontaneous sym- 
metry breaking. The DM candidates are assumed to be 
heavy compared with the vector and Higgs bosons and 
the SM particles are uncharged under this new gauge sec- 
tor. Hence, all the interactions with the SM are mediated 
by kinetic mixing of U(l)s with the photon, if we neglect, 
as usually done, the mixing with the Z boson for the pro- 
cesses of interest here. The Lagrangian can then be written 



C 



_ r-idark rp^v. 
■4 P v 



dark 



_ zTidark rpliv 



\D^\ 2 -V(<P) (2) 



where Fy V is the photon field strength, Fj}„ is the U(l)s 
(dark photon) field strength, and the covariant derivative 
is Dy = d^ + ie'A'^ with 11(1)5 charge e' . The Higgs poten- 
tial V((j)) is assumed to be of a form which spontaneously 
breaks the U(l)s symmetry, with a vev ((/)) = v' /y/2. Af- 
ter spontaneous symmetry breaking, the U(l)s boson ac- 
quires a mass Mjj given by Mjj = e'v' in the unitary 
gauge. As noticed in the literature [26], one can remove 
the kinetic mixing term of Eq. (pQ) by redefining the pho- 
ton A^ — >• Ay — eA^ so that the SM fermions pick up a 
small U(l)5 charge ~ ee, or, alternatively, one can treat 
the kinetic mixing term as an interaction as long as e is 
small. 

With this model at hand, we performed an exact low- 
est order calculation of the radiative processes e + e~ — » 
/i + /i~7, e + e~7, considering all possible s and t-channel 7 
and U boson exchanges and one-photon emission from any 
leptons (see Fig. [I]). Therefore, at a difference with respect 
to earlier works, e.g. Ref. [26], the present calculation in- 
cludes (i) finite width effects for s-channel annihilation 
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subprocesses, both for muon and electron final states and 
(ii) non-resonant t-channel U boson exchange and s-t in- 
terference contributions for electron final states. As far as 
the former ones are concerned, they are of the order of 
F\j jMjj on the integrated cross section and as such fairly 
small for the resonance of interest here. However, the finite 
width effects are a crucial ingredient from a phenomeno- 
logical point of view because they allow, at a variance of 
an on-shell calculation, a direct evaluation of the differ- 
ential cross section of main experimental concern, i.e. the 
invariant mass distribution (see Fig. [2} in the presence of 
arbitrary cuts, as well as an exclusive generation of the 
momenta of the particles coming from the U boson decay. 
Moreover they are crucial in order to reconstruct the cor- 
rect shape of the resonance, in particular the dip induced 
by 7(5) — U(s) interference, as shown in Figs. [3] and HI 
Concerning ^-channel contributions, we checked through 
a number of numerical experiments that their impact is 
rather moderate, being at the per mille level or below it 
for typical e values in the 10 -3 range, in agreement with 
the naive expectation of its order of magnitude. Anyway, 
they can in principle be useful to analyze non-resonant U 
boson production in the ^-channel by studying differen- 
tial distributions such as angular distributions or asym- 
metries in the electron channel, that can be exploited in 
principle as a complementary search strategy, whenever a 
future flavor factory will be able to collect enough statis- 
tics. In our approach the pure U boson signal contribu- 
tion turns out to be defined as the difference between the 
full and background matrix elements. The calculation has 
been done using and adapting (to U boson couplings) the 
algorithm ALPHA [31], which is an efficient tool to com- 
pute tree level matrix elements involving many Feynman 
diagrams. We implemented it in the Monte Carlo genera- 
tor BabaYaga@NLO [32H33ll34lf35l l36] . a standard tool for 
luminosity measurement at flavor factories [30], to sim- 
ulate distributions of experimental interest and account 
for realistic event selection criteria. In the U boson prop- 
agator we include the total dark photon width using the 
formulae of Ref. [25] for the partial widths into leptons 
and hadrons, i.e. 
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>hadrons/«e+e-^/i+/i- is taken into ac- 
count according to the compilation of Ref. [37] and fy = 



1 
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It is worth noting that the U boson is a dra- 



matically tiny resonance, its total width being in the range 
r v ~ 10- 5 - l(T 2 MeV and r v ~ 1(T 7 - 1(T 4 MeV for 
e = 10 -2 and 10 -3 , respectively, and M\j ~ 0.1 - 1 GeV. 
The above feature can be clearly seen in Fig. [2j which 
shows the invariant mass distribution of the muon pairs for 
three values of the kinetic mixing parameter and U boson 
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Fig. 2. Invariant mass distribution of the muon pairs for differ- 
ent values of the kinetic mixing parameter and U boson mass 
at a ^-factory, with y/s = 1.02 GeV. Predictions at the tree 
level. 



mass at DA^NE energies (y/s = 1.02 GeV), as obtained 
through our calculation in the lowest order approxima- 
tion. The contribution of the QED background processes 
is taken into account by means of the exact matrix ele- 
ments for e + e~ — » 7* — » / + /~7, I = e, /iQ implemented in 
the current version of BabaYaga@NLO [35,36 . 

At flavor factories multiple soft and collinear radiation 
emitted by the colliding beams may have a strong impact 
on the measured cross section and on the shape of the dis- 
tributions. The effect of higher order corrections is taken 
into account using the popular QED structure function ap- 
proach [38,39 . Initial state radiation (ISR) modifies the 
tree level cross section as follows 

f 1 , ^ 

da(s,t) = / dxidx2dao(s,t) D(xi,s) D(x2,s) (5) 



where s, t are the Mandelstam invariants after initial state 
photon radiation and the electron structure function D(x, s) 
describes multiple soft and hard photon emission in the 
collinear approximation according to the expression 
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where /3 = 2a/ir (\og(s / m 2 e ) — l) is the large collinear 
factor dependent on the radiating fermion mass and x the 



1 For the electron final state, there is another relevant back- 
ground which should be considered, namely e + e~ — >• 77 with 
conversion of one of the two photons in the beam pipe or the 
inner wall of the detector. This background would require ap- 
propriate studies because of its instrumental origin and can 
be, however, made negligible for small dark photon masses by 
requiring cuts on the reconstructed invariant mass and vertex 
of the pair [28ll2l] . 
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Fig. 3. The invariant mass distribution of muon pairs, for 
standard cuts at DA#NE (y/s = 1.02 GeV), and dark photon 
parameters Mjj — 0.98 GeV, e = 1 X 10 -3 . Results are shown 
at the lowest order (dotted line), including ISR only (dashed 
line) and ISR + FSR (solid line). 



four- momentum fraction of the electron (positron) after 
photon shower. In Eq. (J6j) F is the Euler gamma- function 
and je ~ 0.5772 the Euler-Mascheroni constant. Because 
typical event selections are characterized by cuts on the 
kinematical variables of the observed particles, the in- 
clusive conditions of the Kinoshita-Lee-Nauenberg (KLN) 
theorem for the cancellation of logarithmic enhancements 
are not fulfilled and final state radiation (FSR) must be 
taken into account as well. We implement it in our cal- 
culation by including two further structure functions in 
Eq. ([5]) to describe photon emission from the final state 
leptons, according to the following equation: 

da(s, t) — I dx\ dx2 dy\ dy 2 d(j$(s, i) 
Jo 

D e (x u s) D e (x 2 ,s) D l (y u s') D l (y 2 , s'). (7) 

In the equation above the four structure functions D are 
given by the explicit expression of Eq. (J6j). In particular, 
D e (xi,s) are the structure functions of initial state elec- 
trons giving the probability of finding, inside an electron, 
an electron with four-momentum fraction xi at a virtu- 
ally scale s, D l (yi,s f ) are the structure functions of the 
final-state leptons, computed with (3 = f3f = j3 ei j3^ for the 
electron and muon channel, respectively, and at the scale 
s' given by the lepton pair invariant mass. Equation ([7]) 
takes into account any experimental cuts, that are imple- 
mented by means of a rejection algorithm on an event by 
event basis. More details on the approach can be found in 
Ref. [|0]. 

Being an effect of size comparable to that of photon ra- 
diation, we also consider in our calculation the running of 
the electromagnetic coupling constant in the hard scatter- 
ing cross section according to (see Ref. [30] and references 
therein) 

a(q 2 ) = 
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Fig. 4. The same as Fig. 3 at a B/super-B factory (y/s 
10 GeV), and dark photon parameters Mu — 6 GeV, e 
5 x 10" 3 . 



Aa(q 2 ) = Aai(q 2 ) + Aa h (q 2 ) 



(8) 



1 - Aa(q 2 



where Aai(q 2 ) is the leptonic contribution to photonic 
vacuum polarization (analytically known) and Aah(q 2 ) is 
the non-perturbative hadronic contribution, included ac- 
cording to the parameterization of Ref. [H]. The typical 
effect of vacuum polarization is to enhance the cross sec- 
tion of some per cent, as explicitly checked in the present 
calculation. 

In order to clarify the role of ISR and FSR in the peak 
region, we show in Fig.[3]the impact of ISR and FSR sepa- 
rately on the invariant mass distribution of the muon pair 
in a narrow region around U boson resonant production. 
A U boson mass of 0.98 GeV at y/s = 1.02 GeV is con- 
sidered. Because the Breit-Wigner shape of the muon pair 
invariant mass is induced by the emission of a sufficiently 
hard, detected photon (radiative return mechanism), it 
can be a priori expected that the additional ISR is dom- 
inated around the peak by the emission of soft photons 
for a U boson mass close to the available cm. energy. 
Actually, not to destroy the resonance, ISR photons emit- 
ted additionally to the radiative return photon must have 
an energy fraction smaller than 1 — Mjj/y/s = AE/E, 
AE/E being <C 1 if M v ~ y/s. Therefore, ISR must have 
the typical soft photon effect of reducing the peak cross 
section, as observed e.g. in Z resonant production at LEP. 
This effect is clearly visible in Fig. [3] On the other hand, 
independently of the U boson mass value w.r.t. the cm. 
energy, the main contribution of FSR is degrading the mo- 
menta of the final state leptons. This reduces the fraction 
of events with invariant mass peaked at the U boson mass 
and enhances the number of signal events with invariant 
mass smaller than the resonance value, giving rise to a tail 
on the left of the peak. These effects due to FSR can be 
clearly seen in Fig. [3] 

The separate contribution of ISR and FSR is also shown 
in Fig. 1 for Mu = 6 GeV at y/s = 10 GeV. The above 
reasoning about FSR still applies and a left-tailed broad- 
ening of the peak is again observed. On the contrary, be- 
cause now Mu is definitely smaller than the available cm. 
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energy, no tight cut off is present for the maximum energy 
carried away by the ISR photons. Said differently, resonant 
production is in this case possible also in the presence of 
sufficiently hard IS photons whose effect is, through con- 
volution, enhancing the peak. This can be clearly observed 
in Fig. |U 

More in general, we studied the impact of photon ra- 
diation both for the muon and electron channels, for dif- 
ferent values of the U boson mass and kinetic mixing pa- 
rameter. An example of the effects due to higher order 
QED corrections, as obtained through our generator, is 
given in Fig. [3j which shows the relative contribution of 
the radiation on the invariant mass spectrum of the lepton 
pairs in the presence of standard selection cuts at DA^NE, 
y/s = 1.02 GeV and taking also into account detector res- 
olution effects (see the next section). It can be seen that 
the differential cross section of the muon final state re- 
ceives a moderate positive correction, of some per cent, 
for a relatively small muon-pair invariant mass, (much) 
smaller than the available centre of mass (cm.) energy, 
while a larger negative correction, at the level of 10-20%, 
shows up for M^^- c± y/s. For the electron channel the 
steeper behavior of the lowest order invariant mass dis- 
tribution makes the corrections much larger, of the order 
of +20% for small/intermediate invariant masses and of 
about -30-40% for M/+/- ~ y/s. In the figure a U boson 
mass of 0.98 GeV has been considered, but the shape of 
the correction is substantially independent of the choice 
of cm. energy and Mjy, and thus reflects a general situa- 
tion. If one considers the constraints on the U boson mass 
coming from astrophysical observations, i.e. Mjj < 2M P , 
M p being the proton mass, in order to avoid antiproton 
production through DM annihilation, the latter situation 
is of interest only for searches at the ^-factory DA^NE 
and can be explained in terms of the dominance of soft 
photon emission when Mj+j- ~ y/s. In a broader class of 
new physics models, this relevant effect shows up when 
the U boson mass is close to the kinematical boundary 
imposed by the available cm. energy. As a whole, these 
results point out that the control of photon radiation ef- 
fects is mandatory at flavor factories for a reliable simu- 
lation of all the distributions useful as discovery tools of 
a dark photon signature. 



3 Experimental Sensitivity 



On the grounds of the calculation described in Section [2j 
we revisited the experimental sensitivity to a dark force 
signal evaluated in the literature, and considered novel 
search strategies and event selections not addressed so far. 
For concreteness, we consider the case of the KLOE/KLOE- 
2 experiment at the upgraded DA^NE [42] and of present 
and future experiments at B and super- B factories |43l 
[33]. We assume an integrated luminosity L of 5 fb~ for 
KLOE/KLOE-2 and L = 500 fb"\ 100 ab" 1 for a B and 
super- i? factory, respectively. We compute the statistical 
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Fig. 5. Relative effect of photon radiation on the invariant 
mass distribution (for brevity da in the plot) of the muon and 
electron pairs, for standard cuts at DA^NE (y/s = 1.02 GeV), 
and dark photon parameters Mjj = 0.98 GeV, e — 1 x 10 -3 . 
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requiring the above ratio to be greater than five for dis- 
covery. In Eq. © dp is the full cross section including the 
exchange of virtual photons and U bosons, gb the back- 
ground cross section, Ns and Nb the expected number 
of signal and background events, respectively. To simulate 
detector acceptances, we also impose the following energy 
and angular cuts (from now on large angle selection) 

KLOE 35° < Z±>7 < 145° E l±n > 10 MeV 

B factory 30° < /±>7 < 150° # 7 ,(/±) > 20, (30) MeV 

using as cm. energies y/s = 1.02, 10.56 GeV, respectively. 
Additionally, signal events can be detected as peaks in 
the lepton pair invariant mass close to the value Mjj, as 
visible in Fig. [2j but only in a window Mjj ± 5m, where 
5m is a cut to reject backgrounds or the mass resolution 
of the detector. To avoid washing out the signal against 
the background, 5m should optimally coincide with the 
detector resolution, which is a crucial parameter for these 
searches. As detector resolutions we use 5m = ±1 MeV for 
KLOE/KLOE-2 and the values obtained according to the 
empirical relations of Ref. [26] for B/super-B factories, 
giving 5m ~ =t [1 — 10] MeV for a mass Mjj starting from 
0.1 GeV and up to a few GeV. 

In Fig. [6] we show the reach potential of KLOE/KLOE- 
2 experiment whose sensitivity has been discussed in the 
literature only in the exploratory work of Ref. [23]. We 
give results for the two leptonic final states according to 
our predictions in the lowest order approximation. Ac- 
tually, we observed that the signal and background cross 
sections are both affected by corrections of about the same 
amount; hence the systematics induced by photon radia- 
tion largely cancel in the experimental sensitivity and the 
conclusions in the presence of QED radiative corrections 
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confirm the ones obtained in the lowest order approxi- 
mation. The sensitivity to the kinetic mixing parameter 
e is shown as a function of the U boson mass. The grey 
areas correspond to the exclusion limits imposed by the 
precision measurements of the anomalous magnetic mo- 
ment of the leptons and of q^qed [I? • As already em- 
phasized in previous studies, the muon channel has bet- 
ter reach than the e + e~ channel which is affected by a 
large radiative Bhabha scattering background. Both chan- 
nels have a sensitivity which is significantly degraded if 
the U boson mass is around the p resonance because the 
branching fraction U — >■ / + /~ is suppressed by the dom- 
inant decay mode U — >> 7r + 7r - . Therefore, in the region 
Mu around the mass of the p it would be more conve- 
nient to utilize the 7r + 7r - decay mode to recover the loss 
in the reach due to suppression of the leptonic branching 
ratio. We note that the maximum sensitivity achievable at 
the upgraded DA^NE with a luminosity L ~ 5 fb _1 , i.e. 
e ~ 0.001 — 0.002, is equivalent, as we explicitly checked, 
to that of the present 5-factory experiments BaBar/Belle 
with L ~ 500 fb _1 [26 , since the reach on e follows, from 
Eq. (|9j), the rule e 2 ex (s/L) 1 / 2 . For the case already ad- 
dressed in the literature of a 5-factory with an integrated 
luminosity of 500 fb _1 we compared our results for the 
reach potential with those derived in Ref. [26], observing 
very good agreement. This demonstrates that the sensi- 
tivity limits obtained considering a U boson signal in the 
lowest order, on-shell approximation are robust and un- 
affected by a full signal calculation, since, as argued in 
Ref. [56] , the interference between the resonance and con- 
tinuum background is negligible because of the extremely 
small width of the U boson. Furthermore, the sensitivity 
is not significantly affected by the contribution of QED 
corrections, especially in the muon channel for small or 
intermediate values of the invariant mass, as can be qual- 
itatively understood as follows. Because the tree level and 
QED corrected signal and background cross sections re- 
ceive relative corrections 5 of about the same size, when 
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Fig. 7. The same as [6] at a super-B factory with an assumed 
luminosity L — 100 ab _1 . 



considering realistic detector resolution values that tend 
to reduce the impact of the radiation on the signal cross 
section, the following relations hold 



r^^(l + S)a° s aT D ^(l + S)a% 



(10) 



where, as above, the signal cross section as is defined as 
the difference between the full and background cross sec- 
tion. It then follows from Eq. (J9j) that the maximum sen- 
sitivity achievable according to a lowest order prediction 
is modified by the effect of photon radiation as 



e ^ e QED 



6° 



VTTsj2e° 



(11) 



that explains the results obtained in our simulations. This 
a posteriori reinforces the robustness of the lowest order 
sensitivity results, showing that they are stable against 
QED radiative corrections. Note that, according to Eq. ([TT]) 
and the numerical results shown in Fig. [3j the experimen- 
tal sensitivity is deteriorated by QED corrections when 
multiple soft photon emission dominates, like for Mu ~ 
y/s, whereas it slightly improves, especially in the elec- 
tron channel, for small and intermediate U boson mass 
values, because of the positive correction induced by hard 
bremsstrahlung. 

In Fig. [7] we present the discovery potential at a super- 
B factory, a situation not yet explored in the literature, 
showing again the results for the two lepton final states. 
The muon channel has a definitely better reach than the 
electron one, like at DA^NE. The sensitivity that can be 
reached at a super- B is well below 10 -3 for all mass val- 
ues, with the exception of the region Mu ~ M^. However, 
the latter situation corresponds to a very narrow interval 
because of the extremely small value of the branching ra- 
tios of the (j) resonance into leptons. The same holds true 
for all the other narrow hadronic resonances, not shown 
in the figure because of their tiny width. Hence the large 
data set which can be in principle collected at a super- £? 
collider will allow to probe values of the kinetic mixing pa- 
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Fig. 8. Discovery potential for the /i + /i~7 channel as a func- 
tion of the dark photon mass Mu, at KLOE/KLOE-2 for an 
integrated luminosity L = 5 fb _1 . The sensitivity according 
to the small angle selection cuts as in the text is shown in 
comparison with the large angle selection. 



rameter about an order of magnitude smaller than those 
reachable by present flavor factories q 

In Figs. [8] and [9] the sensitivity for the cuts defined as 
(from now on small angle selection) 

KLOE 35° < 6i± < 145°, E l± > 10 MeV 

|cos(0 7 )| >cos(15°) 

is shown in comparison with the sensitivity for the large 
angle selection previously discussed. This event selection is 
actually used at DA^NE for the measurement of the ttttj 
cross section [30] and has been never considered before 
for this kind of studies. The results are shown for lepton 
pairs invariant masses larger than 0.5 GeV, since below 
this value the cross sections drop steeply to very small 
values. As can be noticed, in both the electron and muon 
channel there is a relevant gain in sensitivity, making this 
selection a very interesting research tool. 



4 Summary 

To summarize, in the present paper the search for a light 
and weakly coupled new gauge boson at contemporary and 
future flavor factories has been analyzed. After perform- 
ing an exact tree level calculation of the signal and back- 
ground processes contributing to the signature e + e~ — )> 
/i + /i~7, e + e~7, we computed the effect of the most im- 
portant higher order corrections induced by multiple pho- 
ton radiation and vacuum polarization, that are a source 



2 We notice in passing that the sensitivity in the electron 
channel could be significantly improved by adopting asym- 
metric angular cuts for the electrons, namely by cutting the 
forward peaking behavior due to t-channel photon exchange, 
thus increasing the signal/background ratio. Nevertheless, the 
muon channel remains the most performing in the large angle 
selection. 
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of relevant systematic effects. Next, we used our calcu- 
lation in order to assess its impact on the experimental 
sensitivity to the dark photon parameters as evaluated in 
the literature, and showed how this can be enhanced by 
means of event selections not considered so far. 

More in detail, we have shown that the effect of QED 
radiative corrections can be of the order of tens of per cent 
on the invariant mass distributions and, more in general, 
non-negligible in the computation of the various observ- 
ables of physical interest. Therefore, in a sensible data 
analysis, these corrections cannot be ignored. 

Using our calculation we have revisited the experimen- 
tal sensitivity to a dark photon signal at present and future 
high-luminosity flavor factories. Moreover we considered 
also new event selections, not explored so far. 

As far as the first item is concerned, we have shown 
that the reach potential is only slightly affected by the 
inclusion of photon radiation effects, and thus a posteriori 
it has been shown to be robust with respect to the inclu- 
sion of potentially dangerous effects. We have shown that 
the maximum sensitivity on the kinetic mixing parameter 
is achievable through the study of the muon channel for 
all values of the U boson mass above the muon pair pro- 
duction threshold, while the electron channel can allow to 
probe the U boson mass region below it with compara- 
ble potential, provided the background e + e~ —> 77, with 
conversion of one photon into an electron pair, is appropri- 
ately rejected. We have analyzed two interesting cases not 
yet addressed in the literature corresponding to the lumi- 
nosity achievable at KLOE/KLOE-2 and a future super- £? 
factory. The statistical significance to the kinetic mixing 
parameter is of the order of e - 0.001 at KLOE/KLOE-2 
and present 5-factories and at the level of a few units in 
10 -4 at a very high luminosity super-!?. 

As far as the second item is concerned, we considered 
a small angle event selection, similar to the one already 
used at DA^NE for the measurement of the ttttj cross 
section, finding that this allows to obtain a considerable 
enhancement in the sensitivity for a U boson mass above, 
say, 0.5 GeV, both for the electron and muon channels. 
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Our calculation is available in an improved version 
of the generator BabaYaga@NLO, which can be used for 
event simulations and data analysis in these new physics 
searches E5l. 
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